Symmetry properties of the large-deviation function of the velocity of a self-propelled 

polar particle 
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A geometrically polar granular rod confined in 2-D geometry, subjected to a sinusoidal vertical 
oscillation, undergoes noisy self-propulsion in a direction determined by its polarity. When sur- 
rounded by a medium of crystalline spherical beads, it displays substantial negative fluctuations in 
its velocity. We find that the large deviation function (LDF) for the normalized velocity is strongly 
non-Gaussian with a kink at zero velocity, and that the antisymmetric part of the LDF is linear, 
resembling the fluctuation relation known for entropy production, even when the velocity distribu- 
tion is clearly non-Gaussian. We extract an analogue of the phase space contraction rate and find 
that it compares well with an independent estimate based on the persistence of forward and reverse 
velocities. 

PACS numbers: 05.40.-a, 05.70.Ln, 45.70.Vn 
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When a particle moves under the influence of a driv- 
ing force through a noisy medium, it occasionally moves 
in the direction opposite to the force. The probability 
of such entropy consumption, relative to production, has 
been shown to obey the well-known fluctuation relations 
[l[ , which are a symmetry property of the large-deviation 
function (LDF) 3 of the entropy production rate. Ex- 
periments on a surprising range of nonequilibrium sys- 
tems find behaviour consistent with the FR in- 
cluding athermal systems [l5| where the noise is a con- 
sequence of the driving.. 

In this paper, we study the velocity statistics of a ge- 
ometrically polar particle in a dense monolayer of beads 
on a vertically agitated horizontal surface. The contin- 
uous input of energy through mechanical vibration is 
balanced by dissipation into the macroscopic number of 
internal degrees of freedom of each particle. Agitation 
feeds energy into the tilting vertical motion of the po- 
lar particle, which transduces it, via frictional contact 
with the base, into horizontal movement in a direction 
determined by its orientation in the plane [HSl. The 
particle thus behaves like a noisy self-propelled object, 
with a statistically significant tendency to move in the 
"reverse" direction, i.e., opposite to its mean direction 
of spontaneous motion. We are interested in large de- 
viations of the velocity. Accordingly, let P{Wr) be the 
probability density of Wr{t) = (1/t) ]l^''[V{t')/ {V)] dt' 
where V{t) = v(f) ■ n(t) with v(<) and n(t) are the par- 
ticle velocity and orientation vector in the plane, and 
denotes an average over the time t. The LDF is 
then F{Wr) = lim^^oo(-l/r) InP(VK^). If Wr were 
the entropy-production rate the FR would read F{Wr) — 

F{-Wr) OC Wr. 

Our experimental results are as follows: (i) The parti- 
cle velocity statistics satisfy a large-deviation principle. 
We are able to extract the LDF, F{Wr), and find that 
the corresponding probability is strongly non-Gaussian, 
with a kink at zero [llj. (ii) The antisymmetric part 
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FIG. I: (a) Geometry of the polar particle: the fat arrow 
indicates the mean direction of "self-propelled" motion, (b) 
A typical experimental screen-shot showing the particle find- 
ing its way through a crystalline medium of 0.8 mm alu- 
minium beads, (c) Comparison of normalized velocity distri- 
butions of polar particle without (gray bars) and with (black 
bars) medium, (d) Mean-square displacement (MSD) with 
and without medium particles. Short-time exponent (with 
medium) and long-time exponents (for both cases) displayed. 



F{Wr) - F{-Wr) OC Wr, i.e., the velocity obeys the 
analogue of a fluctuation relation, (iii) From the veloc- 
ity statistics, we calculate an analogue of phase-space 
contraction rate and show that it correlates very well 
with an independent estimate based on the difference in 
the persistence rates of negative and positive velocities. 
The need of such an independent estimate for contraction 
rates, in a fluctuation-relation context, has been empha- 
sized recently [22 1. 

Two clarifications are essential here: (a) In principle, 
we are not measuring the LDF for the entropy produc- 
tion rate, as the distributions of power and velocity are 
distinct for general time-dependent driving; indeed, we 
have no access to the time-series of the propulsive force. 
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We know of no earlier reports of symmetry relations anal- 
ogous to the FRs for large deviations of the velocity [2^ . 
(b) Such a relation would hold trivially if the velocity 
statistics happened to be Gaussian. We will see that the 
velocity statistics of our particle is far from Gaussian. 

Our experimental cell is a shallow circular geometry, of 
diameter D = 13 cm, made of hardened aluminium alloy. 
To confine the particles to two dimensions, a glass lid is 
fixed on the external perimeter of the circle at a height 
of 1.2 mm above the base. We ensured that base and lid 
were uniformly flat to within 10 ^m accuracy. This ge- 
ometry is mounted on a permanent magnet shaker (LDS 
406/8) and is shaken at a fixed frequency / = 200Hz and 
amplitude ap between 0.019 mm and 0.047 mm. The re- 
sulting accelerations F = ao(27r/)^/(7, measured by an ac- 
celerometer (PCB Piezotronics 352B02) and nondimen- 
sionalized by gravity g, lie between 3.0 and 7.5. We en- 
sured that our apparatus was level to the accuracy of a 
spirit-level; our results are insensitive to small deviations 
in levelling. The dynamics of the particles is recorded 
by a high-speed digital camera (Redlake MotionPro X3) 
mounted vertically above the plate, with a maximum res- 
olution of 1024 X 1280. The frame rate used was 50 fps for 
all the results presented. The resolution of the resulting 
image was 0.1 mm which corresponds to a p ixel length. 
The images were analyzed using Image J 24 1. 

Our "self-propelled" polar particle is a brass rod, 4.5 
mm long and 1.1 mm in diameter at its thick end as 
shown in Fig. [ija). On a bare surface the dynamics 
is as follows: When confined between horizontal plates 
and shaken at F > 4.5 it moves on average with nar- 
row end forward along the arrow in Fig. [UJa), which we 
term the positive velocity direction. We determined the 
position and, thanks to the shape asymmetry of the par- 
ticle, the orientation of the particle in each frame [25| . 
and extracted the instantaneous particle velocity along 
its axis. Gray bars show a typical distribution of nor- 
mahzed velocity {V/{V)) of this particle at F=7.5 in Fig. 
[UJc), clearly showing the tendency towards systematic 
directed motion with exceedingly rare backsteps. Below 
F=4.5 the same particle gradually starts showing a sub- 
stantial number of negative velocity events whose signif- 
icance we will discuss later. 

We then studied the motion of the polar particle when 
the experimental plate is filled with a close-packed mono- 
layer of aluminium beads of diameter d = 0.8 mm. Note 
that the 1.2 mm gap thickness between base and lid leaves 
a clearance of 0.4 mm above the beads, and 0.1 mm above 
the thick end of the polar rod, allowing the play that 
keeps the system "alive". We study the velocity statis- 
tics of the polar particle as it pushes its way through a 
medium at area fraction $ = N{d/D)'^ > 0.8 where N 
is the total number of beads. At these concentrations 
the beads form a triangular lattice (Fig. Hfb)). Black 
bars in Fig. [IJc) show the normalized velocity distribu- 
tion of the polar particle at F = 7.5, $=0.83. Note the 
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FIG. 2: (a) Typical velocity fluctuations of polar particle in 
bead medium, with dashed line indicating mean velocity. All 
data except in (d) is for F = 7.5, ^ = 0.83. (b) Probability 
distribution of Wt for various r. Note exponential negative 
tail and flat peaks, (c) The collapsed large-deviation func- 
tion. Inset shows kurtosis (K) and skewness (S) plotted as a 
function of r clearly indicating the strongly non-Gaussian na- 
ture of the distributions, (d) Large-deviation function shown 
for r = 6.5 and $ — 0.83 and corresponding kurtosis (K) and 
skewness (S). (e) Linear dependence of In P{+Wt) /P{ — Wt) 
on Wt; r=0.04 s, 0.10 s and 0.30 s and (f) Data collapse 
of {l/T)lnP{+Wr)Pi-Wr) vs Wr outo a siuglc line for aU 
T > 0.04s for r = 7.5 and $ = 0.83. 



significant weight at negative velocities. Fig. [Ijd) shows 
the mean square displacement (MSD) plot of the polar 
particle with and without the medium. We see that the 
particle which moved ballistically at all time scales in the 
absence of a medium shows, when surrounded by beads, 
a sub-ballistic short-time motion (MSD ^ t^-^) which we 
attribute to short-time negative-velocity events. 

Typical fluctuations in the velocity of the particle are 
shown in the time-series in Fig. UKa), which has an auto- 
correlation time < 0.02 s. Approximately 50,000 frames 
were captured and the velocity obtained from every pair 
of successive frames. Negative events arc clearly visi- 
ble. After evaluating Wr and corresponding probabil- 
ity P(VFr), our aim is to obtain an LDF F{Wr). To 
this end wc construct the time series of Wr , dividing the 
V{t)/ {V{t)) series into different bins of length r and av- 
eraging over overlapping bins where the centre of each 
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FIG. 3: (a) Typical velocity time series of polar particle with- 
out bead medium, at F = 3.0; dashed line again indicates 
mean velocity, (b) Probability distribution of Wt for r = 
0.06s, 0.22s and 0.70 s. (c) kurtosis (K) and skewness (5*) as 
a function of r indicating tendency to become more Gaussian 
at higher r. (d) {l/r)ln P{+Wt)/ P{-Wr) vs. W^, shows a 
collapse only for r > 0.20s. 



bin is shifted from the previous one by a time difference 
0.02 s to improve statistics. The results presented here 
are not sensitive to the value of the time- difference used 
in the analysis. Fig. [H^b) shows the probability distribu- 
tion for Wr for T 0.04 s and 0.16 s. The distributions 
are highly non- Gaussian with exponential tails on the 
negative side. 

To bring out the non-Gaussian nature, we calculate the 
skewness S = {6W^) jcr" and kurtosis K = i6Wf) ja^ - 
3, where bW^ = Wr - {Wr) and = {5W^). Note 
that S* = for a Gaussian distribution. It is clear 
from the inset to Fig. [5Jc) that the distributions of Wr 
are highly skewed towards positive values at low r, and 
become flatter around the mean at large r. 

In order to extract the LDF from our data, we be- 
gin by writing P{Wr) = exp [— rF(W^i-)], allowing 
for a prefactor Ar independent of Wr and varying more 
slowly than exp(— t). As F is expected 01 to vanish 
at the most probable or typical value of Wt-, we esti- 
mate Ar by the maximum value of P{Wr)- We evaluate 
{-l/T)hY[P{Wr)/Arl find data collapse for r > 0.12s, 
and obtain the LDF. Fig. Hfc) shows F{Wr) for r = 
0.20s, 0.30 s & 0.40 s covering almost the entire range 
of Wr- F{Wr) shows a sharp kink at zero, remaining 
almost flat between and 1. Fig. WiA) shows the LDF 
obtained for yet another case, F = 6.5 at $ = 0.83. The 
distribution is again non-Gaussian though K and S are 
low (see inset in Fig. (Hd)), and a kink at zero can be 
seen. We have no explanation for the difference in the 
behaviours of K in Figs. [l](c) and (d). 

A comparison is in order to the theoretically predicted 



LDF of the entropy production rate for a colloidal par- 
ticle driven by a constant force through a periodic po- 
tential [26[. The resemblance of our LDF (Fig. Wis)) to 
their Fig. 1, bottom row, centre panel, is striking. Note 
however that we are measuring the LDF for the velocity, 
not the entropy production as in [26j . If the motion of 
our polar particle can be approximated as propelled by a 
constant force, then our results can be viewed as a confir- 
mation of the predictions of 26| , albeit in a slightly more 
complicated medium. If it turns out that the propulsive 
force in our case has significant time-dependence, our re- 
sults are all the more intriguing. At present however we 
have no independent way of obtaining a time-series for 
the force on the particle. 

We now examine the relative probabilities of posi- 
tive and negative coarse-grained normalized velocities 
Wr. We find that In [P(-f W^)/P(-W^)] is linear in 
Wr, as shown in Fig. ^e) for t=0.04s, 0.10 and 0.30s 
(only three r values shown for clarity). This linear- 
ity persists to the highest r values where the distribu- 
tion is clearly non-Gaussian. Moreover, Fig. [5Jf) shows 
that {l/T)\TL[P{+Wr)/P{-Wr)] VS. Wr collapscs oiito 
a single straight line for all r values; the slope a ~ 
43.7±5.2 s~^. We find a similar result for other val- 
ues of r and $. We conclude that the antisymmetric 
part of the LDF is linear in W^; i.e F{Wr) - F{-Wr) oc 
Wr- By analogy with the Gallavotti- Cohen SSFR Q 
for entropy flux, we are tempted to suggest a fluc- 
tuation relation for the particle's normalized velocity: 
Vmir^^{^/T)\nP{+Wr)/P{-Wr) = oWr- tt may bc re- 
garded as similar to the phase space contraction rate in 
the conventional SSFR [l[; we return to this point later. 

For comparison we repeat our experiment without the 
bead medium, working at a lower r=3.0 where the parti- 
cle on a bare plate shows frequent negative velocity events 
[13 ■ Fig. [3l[a) shows the velocity time-series and Fig. 
^h) shows the distribution of Wr for r = 0.06s, 0.22s 
and 0.70s. The trend seen here, contrary to that with 
the bead medium, is of a non-Gaussian distribution at 
lower r = 0.06 s, becoming progressively more Gaussian 
as we integrate the time-series for higher r. It is clear 
from Fig. [3Kc) that both skewness and kurtosis decrease 
with r. The quantity (l/rWr) In [Pr{+W)/Pri-W)] ap- 
proaches a r-independent constant only for r > 0.20s 
(see Fig. EKd)) where both non-Gaussian parameters are 
signiflcantly low. This is unlike the case where the par- 
ticle moves through the bead-bed, where this behavior 
persisted even when the distribution was noticeably non- 
Gaussian. 

Returning to the case with a bead-bed, we repeated 
the experiment for various combinations of $ and F. We 
found that in all cases the antisymmetric part of the LDF 
for the velocity was linear, and obtained a range of a val- 
ues. We now suggest that the a can be estimated inde- 
pendently, without reference to the LDF, as follows: We 
extract the probability P+ (t) that a particle moving with 
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FIG. 4: Analogue of phase-space contraction rate a vs. /3_ — 
/3+, the difference in persistence rates of negative and positive 
velocities. Filled squares: $=0.83 (F = 7.5, 6.5, 5.5, 4.5); 
filled circles: $=0.82 (F = 7.5, 5.5); filled triangles: $=0.8f 
(F = 7.5, 6.5, 5.5); filled diamonds: $=0.80 (F = 7.5, 6.5, 
5.5). The values for the case of a bare substrate at F = 3.0 
(open star) and the circular Id channel at F = 7.5 (open circle) 
are also shown. The solid line depicts a = (/?_ —/?+). Inset: 
log- linear plot of P+{t) and P-{t) vs. time t, with slopes /3+ 
and /3_ respectively. 



positive velocity at time continues to do so upto time t, 
and similarly P_ (t) for negative velocity. Each is found 
to decay exponentially, with rates /?+ and /3_ respectively 
as shown in Fig. 3] inset. Since /3_ and /3+ respectively 
measure the mean rates of escape from regions of neg- 
ative (atypical) and positive (typical) velocity, it seems 
plausible that the overall relaxation rate of the system 
should be the difference of the two, i.e. a = /3_ — /3+. 
Indeed, Fig. |4] shows a convincing correlation between a 
and /3- — /3+ for a range of $ and F. The data for a parti- 
cle moving at F = 3.0 on a bare plate, as well as that for 
motion at F = 7.5 in the one-dimensional circular track 
mentioned above, were also analyzed in the same man- 
ner (Fig. S]) and confirm the correlation between a and 
/3_-/3+. 

In conclusion, we have shown that a geometrically po- 
lar particle, when energized by vertical vibration and im- 
mersed in an array of spherical beads, displays frequent 
steps in the sense opposite to its mean direction of spon- 
taneous motion. The resulting velocity distribution is 
highly non-Gaussian, and the large deviation function 
(LDF) shows a kink at zero velocity, as in 26]. Most 
intriguingly, the antisymmetric part of the FDF is linear, 
i.e., the velocity fluctuations obey a symmetry relation 
analogous to those known 1] for the entropy produc- 
tion rate. We provide an independent estimate of the 
analogue of a phase-space contraction rate, in terms of 
persistence rates of positive and negative velocities. 

We thank O. Dauchot for valuable comments. AKS 
and SR respectively acknowledge support from a CSIR 
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